TRANSPORT AND CONCENTRATION PROCESSES IN THE 
MULTIDIMENSIONAL ZERO-PRESSURE GAS DYNAMICS 
MODEL WITH THE ENERGY CONSERVATION LAW 
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Abstract. We introduce integral identities to define (5-shock wave type solu- 
tions for the multidimensional zero-pressure gas dynamics 

pt + V- (pU) = 0, 
{pU)t +W -{pU^U) = 0, 

where where p is the density, U S K" is the velocity, H{x, t) is the internal 
energy, x g M". Using these integral identities, the Rankine-Hugoniot condi- 
tions for (5-shocks are obtained. We derive the balance laws describing mass, 
momentum, and energy transport from the area outside the i5-shock wave front 
onto this front. These processes are going on in such a way that the total 
mass, momentum, and energy are conserved and at the same time mass and 
energy of the moving (5-shock wave front are increasing quantities. In addition, 
the total kinetic energy transfers into the total internal energy. The process of 
propagation of (5-shock waves is also described. These results can be used in 
modeling of mediums which can be treated as a pressureless continuum (dusty 
gases, two- phase flows with solid particles or droplets, granular gases). 



1. Strong singular solutions and pressureless mediums 

1.1. L°°-type solutions. Let us recall some classical results. Consider the Cauchy 
problem for the system of conservation laws in one dimension space: 

r Ut+{F{U))^ = 0, in Mx(0, ex.), 

\ U ^ U°, in M X {t = 0}, 

where F : K™ ^ R"" is called the flux- function associated with HI]); J7° : K ^ R" 
are given vector-functions; U ~ U(x,t) = {ui{x,t), . . . ,Um{x,t)) is the unknown 
function with value in M'", and components Uj{x, t), j = 1, . . . ,m; x E M., t > 0. 

As is well known, even in the case of smooth (and, certainly, in the case of 
discontinuous) initial data U^{x), in general, does not exist any smooth and global 
in time solution of system (|l.ip . As noted in the Evans' book [51 11.1.1.], "the great 
difficulty in this subject is discovering a proper notion of weak solution for the initial 
problem (jl.ip " . "We must devise some way to interpret a less regular function as 
somehow "solving" this initial-value problem" [9l 3. 4.1. a.]. But it is a well known 
that a partial differential equation may not make sense even if U is differentiable. 
"However, observe that if we temporarily assume U is smooth, we can as follows 
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rewrite, so that the resulting expression docs not directly involve the derivatives of 
[/" [HI 3. 4.1. a.]. "The idea is to multiply the partial differential equation in (|l.ip 
by a smooth function Lp and then to integrate by parts, thereby transferring the 
derivatives onto (/?" [H 3.4.1.a.;ll.l.l.]. Following this suggestion we shall derive the 
integral identity which gives the following definition of an L°°- generalized solution of 
the Cauchy problem (|1.1|) : U € x (0,oo);R™) is called a generalized solution 

of the Cauchy problem (jl.ip if the integral identity 

( J (u ■ ipt + F{U) ■ ip^^ dxdt + J U°{x) ■ip{x,0)dx ^0 (1.2) 

holds for all compactly supported smooth test vector- functions ^ : R x [0, oo) — )- R™, 
where • is the scalar product of vectors, and J f{x) dx denotes the improper integral 
J^oo /(^) "This identity, which we derived supposing f7 to be a smooth solution 
makes sense if U is merely bounded" [HI 11.1.1.]. 

Theorem 1.1. (sec, e.g., [S] 11.1.1.]) Let C Mx (0, oo) be a region cut by a smooth 
curve r into a left- and right-hand parts 17 zp. Let us assume that the generalized 
solution U of p.ip is smooth on either side of the curve T along which U has simple 
jump discontinuities. Then the Rankine-Hugoniot condition 

[FiU)]^iyi+[U]^iy2 = 0, (1.3) 

holds along T, where n — (1^1,1^2) is the unit normal to the curve T pointing from 
fl- into 

[F(C/)]'=^F(C/_)-F(C/+), 

def 

[U] = U- - u+ are the jumps in F{U) and in U across the discontinuity curve T , 
respectively. C/zp are respective the left- and right-hand values of U onT . 
IfT = {(a;,t) : x = (^(t)}, where € (7^(0, +00), then 

n=(t/i,t^2)= , \ (1,-0.(0), (1.4) 



and ()1.3|) reads 



where (•) = A(.). 



[F{U)]^ = <P{t)[U]^, (1.5) 



It is well known that HUE L°°(R x (0, 00); R™) is a generalized solution of the 
Cauchy problem (jl.ip compactly supported with respect to x, then the integral of 
the solution on the whole space 

ui.j),.^IWi.)^, ,>o (1.6) 

is independent of time. These integrals can express the conservation laws of quan- 
tities like the total area, mass, momentum, energy, etc. 

1.2. (5-shocks. It is well known that there are "nonclassical" situations where, in 
contrast to Lax's and Glimm's classical results, the Cauchy problem for a system of 
conservation laws either does not possess a weak L°° -solution or possesses it for some 
particular initial data. In order to solve the Cauchy problem in these "nonclassical" 
situations, it is necessary to seek solutions of this Cauchy problem in class of singular 
solutions called S -shocks. Roughly speaking, a S -shock is a solution such that its 
components contain Dirac delta- functions. 

It is customary to assume that a (5-shock wave type solution was first described by 
Korchinski in his unpublished dissertation |13j in 1977. However, in fact, a solution 
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of this type as well as the Rankine-Hugoniot condition for the one-dimensional con- 
tinuity equation were already derived from physical considerations in the book |33[ 
§7, §12] in 1973. Next, in 1979, A. N. Kraiko [M] considered a new type of discontinu- 
ity surface which are to he introduced in certain models of media having no inherent 
pressure and obtained the Rankine-Hugoniot conditions for them. The system under 
consideration in |14j is the zero-pressure gas dynamics described by the system of 
equations: 

(^--).-((^ - »■ 

where p{x,t) > is the density, u{x,t) is the velocity, p{x,t)u{x,t) is corresponding 
momentum, t{x, t) is the internal energy per unit mass, a; g M. The last system can 
be derived from the Euler equations of nonisentropic gas dynamics 

pt + {pu)^ = 0, {pu)t + {pu^ +p)^ = 0, {pE)t + {{pE+p)u)^^Q, (1.8) 

2 

if we set p = 0, where E ~ ^ — h t is total energy per unit mass. 

According to [131 page 502] , to construct a solution for system (|1.7p for arbitrary 
initial data, we need discontinuities which would be different from classical ones and 
carry mass, impulse and energy. As it turned out these nonclassical discontinuities 
are S- shocks. 

The theory of (5-shocks has been intensively developed in the last fifteen years (for 
example, see [2], [4]- [7], [17]- [20], [29]- [31] and the references therein). Moreover, 
recently, in [23], a concept of (5(")-shock wave type solutions was introduced, n = 
1,2,.... It is a new type of singular solution of a system of conservation laws such 
that its components contain delta functions and their derivatives up to n-th order. 
In [24], [27], the theory of S' -shocks was established. The results [24] and [27] show 
that systems of conservation laws can develop solutions not only of the type of 
Dirac measures (as in the case of i5-shocks) but also the type of derivatives of such 
measures. 

The above-mentioned singular solutions do not satisfy the standard integral iden- 
tities of the type ()1.2p . To define them we use special integral identities and derive 
special Rankine-Hugoniot conditions. These solutions are connected with transport 
and concentration processes [1], [5], [23], [23], |28j . 

In the numerous papers cited above (5-shocks were studied for the system of zero- 
pressure gas dynamics: 

pt + V • (p[/) = 0, {pU)t+V ■{pU®U) = 0, (1.9) 

where p = p{x, i) > is the density, U = (mi(x, t), . . . , Un(x, t)) S M" is the velocity, 
V = (gf^, • . • , gf-) J ■ is the scalar product of vectors, (X) is the usual tensor product 
of vectors. 

The system of zero-pressure gas dynamics (|1.9|) has a physical context and is used 
in applications. This system can be considered as a model of the "sticky particle 
dynamics" and was used, e.g., to describe the formation of large-scale structures of 
the universe [25] , [32] , for modeling the formation and evolution of traffic jams [3] , for 
modeling non-classical shallow water flows [8]. Nonlinear equations (in particular, 
zero-pressure gas dynamics) admitting J-shock wave type solutions are appropriate 
for modeling and studying singular problems like movement of multiphase media 
(dusty gases, two-phase flows with solid particles or droplets). The presence of 
particles or droplets may drastically modify flow parameters. Moreover, a large 
number of phenomena that arc absent in pure gas flow is inherent in two-phase flows. 
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Among them there are local accumulation and focusing of particles, inter-particle and 
particle-wall collisions resulting in particle mixing and dispersion, surface erosion due 
to particle impacts, and particle-turbulence interactions which govern the dispersion 
and concentration heterogeneities of inertial particles. The dispersed phase is usually 
treated mathematically as a pressureless continuum. Models of such media were 
discussed in the papers [H]- [16], [21] -[23]. Equations admitting i5-shocks can also 
be used for modeling granular gases. Granular gases are dilute assemblies of hard 
spheres which lose energy at collisions. In such gases local density excesses and local 
pressure falls [10] , [11] . In [10] , [11] , the following hydrodynamics system of granular 
gas 



Pt + (pw)a; = 0, p{ut + uux) = -{pT)^j;, Tt + uT^ = -(7 - 1)Tm2; - ApT^/^, 



was studied, where p is the gas density, u is the velocity, T is the temperature, 7 is 
the adiabatic index, p = pT is the pressure. It was shown that for non-zero pressure 
this system admits a solution which contains a i5-function in the density p. 

1.3. Main results. As it follows from [H]- [H], for modeling media which can be 
considered as having no pressure we must take into account energy transport. In the 
above-cited papers zero-pressure gas dynamics was studied only in the form p.9p . 
Therefore, we need to study (5-shocks in zero-pressure gas dynamics 



where H{x,t) is the internal energy, |J7p = J2k=i This system is obtained by 
adding an energy conservation law to zero-pressure gas dynamics ()1.9p . As distinct 
from (|1.7p it is more convenient for us to consider as a variable H instead of H = pr, 
where r is the internal energy per unit mass. The reason is that since for singular 
solution t{x, t) and p{x, t) must contain (5-functions, it is impossible to define the 
product T{x,t)p{x,t). 

Under the second thermodynamics law it is natural to supplement the system 
(|1.7p with a state equation t = r(T), where T is the temperature. For (jl.lOp the 
natural state equation is H = H{p,T), moreover, H(0,T) = 0. 

In Sec. m we introduce Definition 12.11 of (5-shock wave type solutions for system 
(|1.10p . Next, using this definition, by Theorem 12.11 we derive the corresponding 
Rankine-Hugoniot conditions for (5-shocks (|2.4p . These Rankine-Hugoniot conditions 
are the direct analog of those that were introduced by A. N. Kraiko [l4] . 

In Sec.[3l we show that (S-shocks are related with the transport processes of mass, 
momentum and energy. According to Theorems 13. H [3?2] the mass, momentum and 
energy transport processes between the area outside of the moving 5-shock wave 
front and this front are going on such that the total mass, momentum and energy 
are independent of time. Moreover, the mass and energy concentration processes 
takes place on the (5-shock wave front. 

2. (5-SHOCK TYPE SOLUTIONS AND THE RANKINE-HuGONIOT CONDITIONS 

2.1. (5-shock type solutions. Throughout the paper we shall systematically use 
some results recalled in Appendix 1X1 Let T ~ {{x,t) : S{x, t) = 0} be a hypersurface 
of codimension 1 in the upper half-space {{x,t) : x € M", t € [0, 00)} C R""*"^, S € 
C°°(M" X [0, 00)), with yS{x,t)\g^^ ^ for any fixed t, where V = . . . , 
Let Ft = {x e M" : S{x,t) = O} be a moving surface in M". Denote by v the unit 
space normal to the surface Tt pointing (in the positive direction) from Q.^ ~ {x G 



+ V • (pU) 
{pU)t+y ■ {pU ®U) 








(1.10) 
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S{x, t) < 0} to rj+ = {x G R" : S{x, t) > 0} such that = j = l,...,n. 



The direction of the vector v coincides with the direction in which the function S 
increases, i.e., inward the domain ilf. The time component of the normal vector 
—G = is the velocity of the wave front Tf along the space normal v. 
For system p.lOp we consider the S-shock type initial data 

(C/°(x),pO(a;),77°(x), x G R"; U^{x), xeTo), 

where p°{x) = ^{x) + e°{x)S{To), (2.1) 
H'\x) = i?0(a;) + ;iO(x),S(ro), 

such that C/° G L°°(R";M"), pP,H° G L°°(R";M), e°,h° G C(ro), Tq = {x : 
S^{x) = 0} is the initial position of the (5-shock front, VS'°(x)|^„_p ^ 0, Ug{x), 
X G Fq, is the initial velocity of the (5-shock, (5(Fo) (= S{S^)) is the Dirac delta 
function concentrated on the surface Fq defined by (jA.8[) : 

{S{So), ^{x)) = f ^{x)dTo, V(pG2?(R"), 



dFo is the surface measure on the surface Fq. 

Similarly to [29l Definition 9.1.] we introduce the following definition 

Definition 2.1. A triple of distributions {U,p,H) and a hypersurface F, where 
p{x,t) and H{x,t) have the form of the sum 

p{x, t) = p{x, t) + e{x, t)5{T), H{x, t) = H{x, t) + h{x, t)5{r), 

and U G L°°(M" x (0,oo);R"), p,H G i°°(R" x (0,oo);M), e,/i G C(F), is called 
a d -shock wave type solution of the Cauchy problem (jl.lOp . (j2.ip if the integral 
identities 

p[ip, + U-W^) dxdt + J^e-^== 

fP{x)ip{x,0)dx + I e°(x)^(a;,0)(iFo = 0, 
5ip dT 



pU(^ipt + U ■ Vipj dx dt+ eUs 



St VT+CP 

U"{x)p^{x)ip{x,0)dx + I e°{x)U^{x)ip{x,0)dro = 0, 



^^Hy^{^^^H)U-^^]dxdi 
e|C/5|2 dV 



(2.2) 



St VTTG2 

pO(x)|C/0(x)p , 



hold for all Lp G 2?(M" x [0,oo)). Here / f[x)dx denotes the improper integral 
Jjj„ fix) dx] dT and dTo are the surface measures on the surfaces F and Fq, respec- 
tively; 

C/.-G^-^ (2.3) 
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is the <5-shock velocity, v is the unit space normal to the surface introduced above; 
— G = j^^, 7^ is the 5-derivative with respect to the time variable (jA.5j) : (5(r) is 
the Dirac delta function concentrated on the surface F defined by (jA.SP : 



(<S(5), (^(x,t)) = 



In view of 
derivative: 



ip{x,t) dTt dt ~ j (p{x,t) 
, the (5-derivative in 



dV 



I). 



^/^TG^' 
can be rewritten as the Lagrangian 



5(p dip „dLp dip 



Dw 



2.2. Rankine— Hugoniot conditions. Using Definition 12.11 we derive the 5-shock 
Rankine-Hugoniot conditions for system (|1.10p . 

Theorem 2.1. Let us assume that Q C M" x (0,oo) is a region cut by a smooth 
hypersurface F = {{x,t) : S{x,t) = 0} into left- and right-hand parts fi^ ~ {{^it) ■ 
^S{x^t) > 0}. Let {U,p,H), T be a 5-shock wave type solution of system (|1.10p {in 
the sense of Definition [^?T|) . and suppose that U,p,H are smooth in ft^ and have 
one-sided limits , p^ , on F. Then the Rankine-Hugoniot conditions for the 
(5-shock 



5_ 

St 



SjeUs) 

St 
e\Usl 

2 



Se 
Jt ' 



{eUs ® Us) = 
+ h)Us) = 



{[pU] - [p]Us) ■ V. 
{{pU®U\ ~ {pU\Ui)-v, 



(2.4) 



LL\U 



hold on the discontinuity hypersurface F, where [/([/, p, i/)] = f{U^,p^,H^) — 
f{U^,p'^,H^) is the jump of the function f(U,p,H) across the discontinuity hy- 
persurface F, is the S -derivative (|A.5|) with respect to t, and Vr^ is defined by 

(TOIi . (TOl) . 

Proof. The first two conditions in (|2.4p were proved in [551 Theorem 9.1.]. 

Let us prove the third condition in (j2.4p . For any test function Lp G 15(17) we have 
(p{x,t) = for {x,t) ^ G, G C n. Selecting the test function (p{x,t) with compact 
support in $7^, we deduce from the third identity in (j2.2p that the third relation in 
((rro)) hold in rj±, i.e., 



P\U\' 



H) 



v.(( 



H)U)^0 for {x,t)en'^ 



2 ■ ■ ■ 2 

Now, if the test function ip{x,t) has the support in il, then 



(2.5) 



o-nG 



p\U\' 



n+nG 



+ HjU -Vip] dxdt 
2 



HjU -S/ip] dx dt 
, (P\U? 



H]U -Vip] dxdt 
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Using the intcgrating-by-parts formula, we obtain 



n±nG 



P\U\' 



n±nG 



H] +V- 



P\U\' 



+ H]U] ]ip{x,t)dxdt 



ma 



|V(..t)5|j 



ip{x, t) dT 



n±nGnM" 

where dT is the surface measure on F. Next, adding the latter relations and taking 
into account (|2.5p . we have 



2 ^ (^)J Vl^;, 0) ax 



P\U\' 



H 



G 



{^ + H)u]..)^ix,t) 



dT 



(2.6) 



Next, applying the integrating-by-parts formula (|A.10[) to the second summand in 
third identity (\2.2^ . one can see that 



e\Us\^ , ,\Sip dT 



h 



5t Vl + G^ Jro 



S^fe\Us\^ 
6t 



+ h}'{x)jip{x,Q)dTo 
dT 



where the adjoint operator jr is defined in (jA.lip . Thus 



C'\Us? AScp dT 



h 



6t vrTG2 



S fe\Us\^ 



6t\ 2 



/i + V 



Pt 



e\Us\' 



+ h''{x))ipix,0)dTo 
dT 



VT+G^ 



(2.7) 



Adding (|2.6p and (j2.7p and taking into account p.2p . (|2.3p . we derive 









- 2 











• + 



- Vr 



2 



dT 



0, 



□ 



for ah ip G 2?(fi). Thus, the third relation in ((^^ holds. 

The right-hand sides of the equations in (|2.4p are called the Rankine-Hugoniot 
deficits in p, pi7, and ^^^^^ — h i?, respectively. 

Let a{x, t) be a smooth function defined only on the surface T — {{x,t) : S{x, t) = 
0} which is the restriction of some smooth function defined in a neighborhood of T 
in R" X M. It is easy to prove that 



Vr, • {aUs) = -2/CGa, 



(2.8) 
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where /C is the mean curvature of the surface Tf (see (jA.7|) ). Indeed, according to 

(EU, dMI, m, we have Vr, • (aUs) = ELi = ELi ^^fe + 

^'^Sl-=i 1^7 = — 2/CGa. Here the obvious relation "^ax""* = ^ was taken 

into account. 

Due to (|2.8p . the Rankine-Hugoniot conditions p.4p can also be rewritten as 

^-2/CGe = ([pC/]-[p]C/5) -i., 

- 2ICGeUs = {[pU(g)U]~ [pUps) ■ ly, 
S f^^lUsW A ^^^fe\UsW ,\ (\(P\U? (2-9) 



5t 



H]U 



p\U 



|2 



H 



Us 



Remark 2.1. The Rankine-Hugoniot conditions (|2.4p constitute a system of second- 
order PDEs. According to this fact, for system (|1.10p we use the initial data (|2.ip 
which contain the initial velocity Ug{x) of a 5-shock. This is similar to the fact 
that in the measure-valued solution approach [3], |17| . |18| . |31| the velocity U is 
determined on the discontinuity surface. 

In the direction v the characteristic equation of system p.lOp has repeated eigen- 
values X ~ U ■ So, we assume that for the initial data (j2.1l) the geometric entropy 
condition holds: 

U"+{x) ■ < U^x) ■ < U"-ix) ■ (2.10) 

where z/" = is the unit space normal of Fq, oriented from D,'^ = {.x G M" : 

S°{x) < 0} to n+ ^{xe R" : S°{x) > 0}. Similarly, we assume that for a solution 
of the Cauchy problem (jl.lOp . (|2.ip the geometric entropy condition holds: 

U+{x,t) ■ zy|j,^ < Us{x,t) ■ < U-{x,t) ■ z^Ij,^, (2.11) 

where Us is the velocity (|2.3p of the (5-shock front Tf , is the velocity behind the 
(5-shock wave front and ahead of it, respectively. Condition (|2.1ip implies that all 
characteristics on both sides of the discontinuity Ft must overlap. For t = the 
condition (|2.1ip coincides with p.lOp . 



3. (5-SHOCK MASS, MOMENTUM AND ENERGY TRANSPORT RELATIONS 

The classical conservation laws (|1.6p do not make sense for a i5-shock wave type so- 
lution. "Generalized" analogs of conservation laws (|1.6p were derived in [5], [H], [H] 
for the one-dimensional case, and in |29] for the multidimensional case. Now we de- 
rive these transport conservation laws for the case of system (jl.lOp . 

Let us assume that a moving surface Ft ~ {x : S{x, t) = O} permanently separates 
into two parts flf ^{xeW : ±S{x, t) > 0}, and f^J {x e R" : ±S°{x) > 0}. 
Let ([/, p, H) be compactly supported with respect to x. Denote by 



and 



M{t)^ / p{x,t)dx, m{t)^ / e{x,t)dTt, (3.1) 

Jn:run+ J Ft 



Pit) ^ I p{x,t)U{x,t)dx, p{t)^ / e{x,t)Us{x,t)dTt, (3.2) 
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masses and momenta of the volume fij U il^ and the moving (5-shock wave front F^, 
respectively, dTt being the surface measure on Tf. Let 

Wkinit) = / dx, Wkinit) ^ / dVt, (3.3) 



and 

W,nt{t)= [ H{x,t)dx, W^nt{t)= [ h{x,t)dTt, (3.4) 

be the kinetic and internal energies of the volume D,^ U Vlf and the moving wave 
front Ft, respectively. Here Wkin{t) + Wkin{t) and Wi„t(t) + Wint{t) are the total 
kinetic and internal energies, respectively; Wkin{t) + Wkinif) + Wint{t) + Wint{t) is 
the total energy. 

Theorem 3.1. Let {U,p,H) together with a discontinuity hypersurface F = {{x,t) : 
S{x,t) =0} be a 5-shock wave type solution [in the sense of Definitional^ of the 
Cauchy problem (|1.10p . (j2.ip . where 

p{x, t) = p{x, t) + e[x, t)5[T), H{x, t) = H{x, t) + h{x, t)S{T). 

Let this solution satisfy the entropy condition \2.1\\ . Suppose that [U,p,H) is com- 
pactly supported with respect to x, smooth in fl^ = {{x,t) : ztS{x,t) > 0} and has 
one-sided limits , p^ , on F. Then the following mass and momentum balance 
relations hold: 

M{t) = ~m{t), m(i)>0, Pit) = -p{t): .3 5^ 

M(t) + m(t) = Af(0) + m(0), Pit) + p{t) = P{0) + piO). 

In fact, the proof of Theorem 13.11 coincides with the proof of [551 Theorem 9.2.]. 
The proof of p9l Theorem 9.2.], and, consequently, the proof of Theorem 13.11 are 
based on the volume and surface transport Theorems lA.ll IA.2I and use the first two 
relations in ([27 



Theorem 3.2. Let {U,p,H) together with a discontinuity hypersurface T ~ {{x,t) : 
S{x,t) = 0} satisfy the same conditions as in Theorem 13.11 Then the following 
energy balance relations hold: 

Wk^n{t) + W^nt{t) > 0, Wk^n{t) < 0, ty„t(t) + W^t (t) > 0, W'„t(t) < 0. (3.6) 

Moreover, 

Wkin{t) + Wkinit) = -{Wintit) + mmt(t)), 
Wku^it) + Wkinit) + Wintit) + W,nt{t) = (3-7) 
= Wk,niO) + Wk^niO) + W.ntiO) + W.ntiO)- 

Proof. 1. Let us assume that the supports of Uix,t) and pix,t) with respect to x 
belong to a compact X € M" bounded by OK. Let = H K. By v we 
denote, as before, the space normal to Ft pointing from fl^ to Qf. Differentiating 
Wkin (t) + Wint (t) and using the volume transport Theorem lA.l) we obtain 

Wk.nit) + mmit)= I |(^i^l»^ + i/(x,^))d. 

J K7UK+ Ot \ Z J 



ik;uki 
'pix,t)\Uix,t)\^ 



+ Hix,t)^Vix,t)-i>drt, (3.8) 



where i> is the outward unit space normal to the surface dK^ and Vix,t) is the 



velocity of the point x in K^. 
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Next, taking into account that for x G system (|1.10p has a smooth solution 
(f/±,p±,i/±), i.e., 



0, 



and C/^, p^, are equal to zero on the hypersurface dK^ except Ft, applying 



Gauss's divergence theorem to relation p.Sp . we transform it to the form 

'p-\u-? 



Ft 



p+\U+\' 



H+]U+] dx 



H-]U-j dx 



H 



Us ■ V dTt 



p-\u-\' 

2 
2 

P\U\' , 



H-)U- -lydTt 



H]U 



H+]U+ -v dTt 
P\U\' 



P\U\' 



H 



Us ■ i^dTt 



H 



Us ■ I' dTt, 



(3.9) 



where Us = is the velocity (j2.3p of the (5-shock front Tt- Using the third 

Rankine-Hugoniot condition p.4p . relation p.9p can be rewritten as 



Wk^n{t)+W,ntit) 



Ft 



-( 



e\Us\' 



h]Us] dT 



(3.10) 



Applying the surface transport Theorem I A. 2 1 to the second relations in (|3.3p . (|3.4p 
one can see that the right-hand side of (|3.10p coincides with —Wkin{t) ~ Wint{t)- 
Thus relations ^J^i hold. 

Since > 0, > and the solution {U,p,H) of the Cauchy problem (jl.lOp . 
(|2.ip satisfies the entropy condition (|2.1ip . we have 

{[p\U\'U]~[p\U\']Us)-^ 



{p-\U-\\U- - Us) ■ iy + p+\U+\^iUs - U+) -ly)]^ > 0; 



(3.11) 



{[HU] - [Hps) ■ V = {B-{U- -Us)-v + H+{Us - U+) ■ > 0. (3.12) 

Formulas dSH), ([3ll|) . ((3l^ imply that W^fcj„(t) + PFi„t(t) < 0, i.e., due to dSH]) the 
first inequality in p.6p holds. 

2. In fact, the second inequality in p.6p was proved in pS] . 

Let us calculate w(<). Taking into account formula (|2.8p . due to the surface 
transport Theorem lA. 21 we obtain 



Wkrnit) = (^■^{e{x,t)\Usix,t)\^) +Vt, ' (e(a;, t) 1 1/^ (x, t ) | ^ C/5) 



(^-(e(a;,i)|f/5(^,0l') - 2ICGeix,t)\Usix,t)\^j dTt 



E 



St 



uske- 



5usk 
St 



-2ICGe{x,t)\Us{x,t)f^dTt. (3.13) 
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According to ((^^ and (g^, wc have 

77 usk + - 2/CGe usk = [pukU ■ v] - [puk]Us ■ v, 

St ^Jt (3.14) 

— usk - 2K.Geusk = [pU ■ v] ugk ~ [p]Us ■ vusk, 

where usk^x, t) is the fc-th component of the vector Us, k ~ I, . . . ,n. Now, subtract- 
ing one equation from the other in (j3.14p . we obtain 

e-^ = [pukU ■ v] - [puk]Us -ly -[pU ■ iy]usk + [p]Us ■ vusk- (3.15) 
Substituting equations p.isp into p.l3p . one can easily calculate 
If " 

Wkin{t) =2 (^X! i^P'^kU ■ v] - [puk]Us ■ v)usk 

k=l 

-[pU ■ iy]\Us\^ + [p]\Us\''Us ■ ly) dVt. 

Taking into account that Us ~ Gv, G = — j^^, i.e., usk — Gi^k, k = 1, . . . ,n, we 
rewrite the above relation as 

wk^n{t)^\j {2[p{U ■vf]G-^[pU ■u]G'' + [p]G'')dTt. (3.16) 
Using p.9p and (|3.16p . we obtain 

Wk^n{t) + Wku,{i) = -\ I {[P\U?U ■ V] - [p|C/p]C/5 • U 



- 2[p{U ■ yf]G + i[pU ■ iy]G^ - [p]G^) dTf (3.17) 
The gas velocity U\rt on the wave front Tt is the sum of the normal component U ■ v 
and the component Utan tangential to the surface Tt- Since |J7p|rt ~ {U ■ + Utam 
and G = Us ■ V, one can represent the integrand in p.l7p as 

{m^U ■ v\ - \p\Uf\Us ■ V - 2[p(C/ • vf\G + 'i\pU ■ v]a' - [p]G^) 

= p-{U^,,y{U- ■v-Us-v)+ p+{U+,y{Us -y-U^-v) 

+ p-{U~ -v ^Us- vf + p+{Us -v-U- ■ vf. (3.18) 
Since a solution {U,p,H) of the Cauchy problem (|1.10l) . (|2.ip satisfies the entropy 
condition p. lip and > 0, wc deduce that the expression ()3.18p is non-negative. 
Formulas ([3Tf| . ([3?T8| imply that Wktn{t) + Wktn{t) < 0. Due to ([SJ]) . we conclude 
that the third inequality in (|3.6p holds. 

3. Since U,p,H are smooth in = {{x,t) : ztS{x,t) > 0}, it easy to see that 
for (x, t) S Kf^ the first and second equations in (jl.lOp imply that 

i=i ^ 

Multiplying the both sides of the above equation by and summarizing over k = 
1, 2, . . . , n, we obtain 

(|[/±|2)^+^,±^i^ = 0, {x,t)eKt (3.19) 



3 



According to p.l9p and the first equation in p.lOp 



j^ + La^i^^-. j=0' (x-,t)Gi^, . (3.20) 
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In the end, from (|3.20p and the third equation in (|1.10p we obtain that 

(i/±)^ + V-(if±C/±) =0, {x,t)eK^. (3.21) 

Next, as before, differentiating Wint{t), using p.21|) and applying the volume 
transport Theorem I A. ![ we obtain 

mntit)= [ ^^^^dx+ [ H{x,t)V{x,t)-vdTt, 

Jk-uK+ JdKtUdK+ 

= - [ \7-{HU)dx+[ H{x,t)V{x,t) ■ DdTt, (3.22) 

Jk-UK+ JdK-UdK+ 

where D is the outward unit space normal to the surface dK^ and V{x,t) is the 
velocity of the point x in . 

Taking into account that U^, , are equal to zero on the hypersurface dKf^ 
except Tt and applying Gauss's divergence theorem to p.22p . we transform this 
relation to the form 

W,,,t{t)^'- [ \/ ■ {H-U-)dx- [ V-{H+U+)dx+[ [H]Us-iydrt 
Jk- Jk+ JTt 

H-U^ -v dTt + 1 H+U+ -ly dTt+ / [H] Us ■ v dVt 



= - [ {[HU] - [H]Us) -lydTt, (3.23) 

JVt 

where [/^ = V |p is the velocity (|2.3p of the (5-shock front Ft, = i>|p is the space 

normal to Ft pointing from K^T to K^^ . In view of the entropy condition ()2.1ip 
the inequality ()3.12p holds, and consequently, p.23p implies the fourth inequality in 
(EH). □ 

Corollary 3.1. According to Theorems 13.11 the mass, momentum and energy 
transport processes between the volume outside of the S-shock wave front fl^ U 
ilf = {x £ M" : S{x,t) ^ 0} and the moving S-shock wave front Tt are go- 
ing on such that the total mass M{t) + m(t), momentum P{t) -\- p{t) and energy 
Wkin{t) + Wkin{t) + Wint{t) + Wintif) are independent of time. More precisely the 
mass and energy concentration processes on the moving 5 -shock wave front Tt are 
going on. In addition, the total kinetic energy Wkin{t) + Wkin(t) transfers into the 
total internal energy Wint{i) + Wint{t). 

The inequality W{t) < in (|3.6p reflects the well-known fact that the evolution 
of a solution with shocks is connected with decreasing of the kinetic energy. 

Remark 3.1. Let us suppose that in a finite time t the whole initial mass M(0) 
and energy Wfci„(0) -I- Wi„t(0) may be concentrated on the (5-shock front Tt- Then, 
according to The Rankine-Hugoniot conditions, for t > t, instead of the whole system 
of zero-pressure gas dynamics (|1.10l) we obtain exactly a "surface" version of this 
system 

Vr-(et/i)=0, 



Se 

a 

5{eUs) 



Vr • [eUs ® Us) = 0, 



^(^-)--((^-^)-)-. 

where Us is the velocity of the (5-shock front Tt, e is the surface density of the front 
mass, h is the surface density of the front internal energy. This system is an analog 
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of the initial system of zero-pressure gas dynamies p.lO[) on the (n — 1)- dimensional 
surface Fj. This (n — 1) - dimensional analog also has the same type as the initial 
system, therefore its solution can develop singularities within a finite time, and all 
mass concentrates on the manifold of dimension n — 2, and so on. Thus, it may 
happen that after the finite number of bifurcations the whole initial mass will be 
concentrated at the singular point. 

4. Example of an one dimensional concentration process 

In the ID case we construct an explicit example of the concentration process 
based on another method. Namely, let us consider the data that do not imply the 
(5-shock initially: 

t/0(.T) = ([/--[c/]e(x))x/(x)), 

p\x) = {p- -[p]Q{x))xi{x), (4.1) 
H\x) = {H- -[H]Q{x))xi{x)). 
where U~ , p~ , H~ ,U~^ , p^ , are constants, p^ , p^ , , > 0, \U] > 0, Xii^) 
is the characteristic function of the segment / = [—L,L],L^ 1. Let us note that we 
can apply the standard mollification procedure to obtain functions smooth at the 
points ±L, but here do not need to do it. 

We obtain the solution to the Cauchy problem by means of the free particles 
method [1]: first we assume that the particles do not feel one others and form the 
overlapping domain. Then we switch to the sticky particles model and change this 
overlapping domain to a point where the mass accumulates according to the conser- 
vation of mass and momentum. Now we have to consider the additional law of con- 
servation of energy. Thus, according to ([1]), the free-particles solution {pfp,Ufp) 
to the two first equations to the zero pressure model has the form 

-L + {U- - [U])t <x<{U- - [U])t, 
(U- -[U])t<x< U-t, 
L + U-t > x> U't, 
otherwise, 



PFp{t,x) = 




{U-, -L + (U- -[U])t<x <{U- -[U])t, 

f/- + ^r-^[C/], [U- -[U])t<x<U-t, 
U--[U], L + U-t> x>U-t. 

Outside of the segment [~L + {U^ — [U])t,L + U~t] the solution Upp contains a 
rarefaction wave, however this part of solution does not contribute to the energy, 
since for the domain of rarefaction p = H = Q. The respective solution {p, U) to the 
sticky particles model is 

p{x, t)= p- - [p\Q{x - Xj (t)) + e{t)S{x - Xj (t)), 

U{x,t) = U- ~ [U]e{x-Xj{t)), 

where the position of the singularity Xj (t) is the following: 



[Up] ^i-[Up]2 + [p][U2p^) 

Xj{t)^ j^j t, if [p]^0, (4.2) 

and 

^ ?Rl-m, . E1±E1, „ H^o. («) 

The amplitude of the (5-shock reads as 

e(t) = [Up]t - [p]x,it). 
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The solution induces the foUowing balance of energy. 

For the sake of simplicity we dwell on the latter case. Thus, 

W^nt{t) = {xj{t) + L- U-t)H- + {L+ {U- - [U])t - Xj{t))H+ 

, , H+ + H- , , 
= Wi„t(0) [U]t, 

it)^l{ix,it)+L-U-t)p-iU-r 

+ {L + iU- -[U])t~x,it))p+iU+r) 

Wfcm = ^[U\p I I t, 

Wint = Wint{0) + Wfc,„(0) - Wi„t(t) - Wktnit) - Wfej„(t) = 

We see that Wkm (t) and W^t {t) decrease with a constant velocity, and vanish within 
a finite time, Wkmt increases unless U~ + ^ 0, (in the latter case w^init) = 
0), Wintit) increases in any case. Since we associate the internal energy with a 
temperature, it signifies that the concentration process always entails the heating of 
point of the mass accumulation and cooling-down of the environment to the " absolute 
zero" that relates to the zero internal energy. 

Appendix A. Some auxiliary facts 

A.l. Moving surfaces of discontinuity. Let us present some results from [T^ 
5.2.] concerning moving surfaces. Let Ft be a smooth moving surface of codimen- 
sion 1 in the space M". Such a surface can be represented locally either in the form 
Tt = {x £ M" : S{x,t) = O}, or in terms of the curvilinear Gaussian coordinates 
s — (si, . . . , s„_i) on the surface: 

Xj = Xj{si, . . . ,Sn^l,t), s € M""-^. 

We also consider the surface F — {{x,t) e M"+^ : S{x,t) = O} as a submanifold of 
the space-time K" x R. We shall assume that \7S{x,t)\^^ ^ for all fixed values 
of t, where V = (g§^, • • • , gf-)- Let v be the unit space normal to the surface Ft 
pointing in the positive direction such that = |V5'|i^j, j = 1, . . . ,n. 

Let f{x,t) be a function defined on the surface Ft for some time interval, and 
denote by ^ the derivative with respect to time t as it would be computed by 
an observer moving with the surface. This derivative has the following geometrical 
interpretation. Let Mq be a point on the surface at the time t = t^. Construct 
the normal line to the surface at Mq. At the time t = to + At, At is suSiciently 
small, this normal meets the surface Ff^+At at the point M = M{to + At). Then 
the (5-derivative is defined as 

Sf{Mo,to) ^ f{M)-f(Mo) 

St At^o At ' ^ ' ' 

If As is the distance between Mq and M, then 

G= lim 4^ (A.2) 

At^O At ^ ' 
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is the normal velocity of the moving surface Tf and 

^ = Ihn ^ = Hm = Gz.,, j = l,...,n. (A.3) 

St At^a At At~^Q At As j> J > > \ J 

Since it is essential that the (5-derivative is computed on a surface, and S remains 
constant on this surface, then ^f- = 0. Thus we have 

SS dS 5S Sxj dS v^^,„„, 9 

j=i j=i 

i.e., 

St^-G\\/S\. (A.4) 

From this formula we can see that — G = can be interpreted as the time 

component of the normal vector. 

The space-time unit normal to the surface F is given by n = ^i+e^ ' 

VTTg^=^^, V(.,) = (v,^). 

If f{x, t) is a function defined only on F, its first order J-derivatives with respect to 
the time and space variables are defined by the following formulas [l2l 5. 2. (15), (16)]: 

Sldg dj^ ^dj ^ _ ^ dj . ^ 

6t dt dv ' 5xj dxj ■' dv' 

where / is a smooth extension of / to a neighborhood of F in R" x M, j = 1, . . . , n, 

and ^ = ^ ' V/ is the normal derivative. Since (5-derivatives are independent of the 
way of extension of the function /, we shall drop tilde from the function /. Thus 
the gradient tangent to the surface Ft is defined as 

Vr,=V-V.= (-^,...,-^), (A.6) 

where 'SI v = v{y ■ V) is the gradient along the normal direction to Ft. The mean 
curvature of the surface Ft is defined as 

i=i ■' 

A. 2. Distributions defined on a surface. The Heaviside function iJ(S') is intro- 
duced by the following definition: 

{H{S), ip{-,-)) = f ipix,t)dxdt, V¥>g2?(]R" xR). 

According to [T2j 5. 3.(1), (2)], we introduce the delta function 6{S) on the surface F: 

{6{S), ^(•,-))= r I ipix,t)drtdt= [ ^{x,t)^=£=, (A.; 



for all Lp e P(R" X R), where c?Ft and dT arc the surface measures on the surfaces 
Ft an F, respectively. According to [TH 5. 5. Theorem 1.], we have 



16 



S. ALBEVERIO, O. S. ROZANOVA, AND V. M. SHELKOVICH 



A. 3. An integration-by-parts formula. Wc need the following iiitcgrating-by- 
parts formula. 

Lemma A.l. ( [291 Lemma 9.1.], cf. [IS 5. 2. (25), (26)]) Suppose that a{x,t) is a 
smooth function defined only on the surface F = |(x,t) : S{x,t) = 0} which is the 
restriction of some smooth function defined in a neighborhood of F in M" x R, and 
Fq = jj: : S{x, 0) = O}. Then the following formula for integration by parts holds: 

f 6ip dT f S*a dT f , . , . . 

for any ip G I?(R" x [0,cx))), where is the adjoint operator defined as 

^^|-2^Ga=| + Vr..(aG.), (A.ll) 
K- is the mean curvature (jA.7|) of the surface Ft . 

A. 4. Transport theorems. Here we give the following transport theorems. 

Theorem A.l. ( [12i 12.8.(3)]) Let f{x,t) be a sufficiently smooth function defined 
in a moving solid fit, and let a moving hypersurface dflt be its boundary. Let v be 
the outward unit space normal to the surface dQt and V{x,t) be the velocity of the 
point X in fit. Then the volume transport theorem holds: 

4 / f{x,t)dx= [ %dx+ [ fV-vdTt 
where dTt is the surface measure on the moving surface dflt- 

Theorem A. 2. ( [T2j 12.8.(9)]) If e{x,t) is a smooth function defined only on the 
moving surface Tt = {x : S{x,t) = 0} {which is the restriction of some smooth 
function defined in a neighborhood ofTt), then the surface transport theorem holds: 

d f 

-j^e{x,t)dTt 

where Ug = vG is the velocity ofTt given by (|2.3p . 
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